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Recently, we have numerically evaluated the current-current correlation function for the ground 
states of three-band models for the CuO planes of high-T c superconductors at hole doping x = 1/8 
using systems with 24 sites and periodic boundary conditions. In this article, the numerical analysis 
is explicated in detail and extended to a wider range of parameters. Our results show no evidence for 
the time-reversal symmetry violating current patterns recently proposed by Varma. If such current 
patterns exist, our results indicate that the energy associated with the loop currents must be smaller 
than 5 meV per link even if the on-site chemical potential on the oxygen sites, which is generally 
assumed to be of the order of 3.6 eV, is taken to zero, as advocated by Varma. We also vary the 
inter-atomic Coulomb repulsion scale and find only a weak dependence on this parameter. So while 
our studies do not rule out the existence of such current patterns, they do rule out that quantum 
critical fluctuations of these patterns are responsible for phenomena occurring at significantly higher 
energies such as the superconductivity or the anomalous properties observed in the strange metal 
phase provided the CuO superconductors are adequately described by any of the three-band models 
discussed. 

PACS numbers: 74.20.Mn, 74.72-h, 74.20.-z 



I. INTRODUCTION 

High-r c superconductivity (HTSC) has been one of the 
most active fields of research in condensed matter physics 
in the past two decade*^. It has turned out to be an 
exceedingly difficult problem, with much of the effort in- 
vested just deepening the mysteries, but it has also led 
to a plethora of new developments extending far beyond 
the field. Many ideas, even though too general to qual- 
ify as complete theories of the cupratcs, have inspired a 
vast amount of research in both high-T c and other ar- 
eas. Most prominently among them are the notions of a 
resonating valence bond (RVB) stated, the gauge theo- 
ries of antiferromagnetismi, and the notion of quantum 
criticality^. There have been, however, a few concise 
proposals which make falsifiablc predictions. Intellectual 
masterpieces among them have been the theory of anyon 
superconductivity^, the proposal of kinetic energy sav- 
ings through interlayer tunneling^, the SO (5) theory of a 
common order parameter for superconductivity and mag- 
netism^, and a more recent proposal that the anomalous 
properties of the cuprates may be due to quantum criti- 
cal fluctuations of current patterns formed spontaneously 
in the CuO planes*^. In a recent Letter—, we investi- 
gated this proposal by finite cluster calculations. Here, 
we provide supplemental information and a more elabo- 
rate account of the approach. 

This paper is organized as follows. In Section [Til we 
briefly review Varma's proposal of spontaneous T vio- 
lation in the CuO planes and discuss some assumptions 
made therein. In Section llTTl we introduce the three-band 
model Hamiltonian we investigate. In Section IIV1 we 
compute the current-current correlations of the ground 
state which we use to obtain information about the ex- 
istence of orbital currents and the magnetic moment as- 



sociated with them. We find that if an orbital current 
phase exists in the cuprates, the energy associated with 
the spontaneous currents will not be sufficiently high for 
the phase to account for the strange metal phenomenol- 
ogy in the cuprates. We derive an upper bound for the 
magnetic moment per unit cell from the upper bound 
we obtain for spontaneous currents, and find it smaller 
than the magnetic moment measured in a recent neu- 
tron scattering experiment^. The comparison shows, 
however, that even if the observed magnetic moments 
were due a current pattern as proposed by Varm a 9 : 10 , 
the magnitude of these currents would be insufficient to 
determine the phase diagram. In search for an alterna- 
tive explanation of the experiment, we investigate the 
spin-spin correlations of the ground state in Section [V] 
In Section IVH the results for various three-band model 
parameters are presented. In particular, we vary the on- 
site chemical potential on the oxygens e p from previously 
3.6 to 1.8, 0.9, 0.4, and finally eV. Furthermore, we 
vary the Coulomb interaction scale from 1.2 to 2.4 
eV. The conclusion regarding the relevance of an orbital 
current pattern for the strange metal phase of CuO su- 
perconductors we reached previously remain intact. In 
Section l"VIH we discuss the the role of finite size effects 
in our numerical experiments, with particular emphasis 
on the net spin 1/2 of our finite size ground states. In 
Section [Villi we conclude that while we cannot rule out 
that the orbital current phase exists in the cuprates, we 
can infer that the energy associated with these fluctua- 
tions is not sufficiently high to account for the strange 
metal phase in the cuprates. 



II. HYPOTHESIS OF SPONTANEOUS T 
VIOLATION IN THE CUPRATES 

The proposal of a spontaneous symmetry breaking 
through orbital currents is motivated by experiment. The 
normal state of the cuprates at optimal doping shows a 
behavior which can be classified as quantum critical, and 
has been rather adequately described by a phcnomcno- 
logical theory called marginal Fermi liquicU 3 -. The linear 
temperature dependence of the normal-state resistivity in 
optimally doped LSCO and YBCQi^±&, which persists 
over several magnitudes of temperature, provides strik- 
ing evidence in favor of this picture. The marginal Fermi 
liquid phenomenology led Varma to assume a quantum 
critical point (QCP) at a hole doping level of a; c w 0.19, 
an assumption consistent with a significant body of ex- 
perimental dat a 16 ' 17 i 18 ' 19 i 20 . Critical fluctuations around 
this point are then held responsible for the anomalous 
properties of the strange metal phase, and provide the 
pairing force responsible for the superconducting phase 
which hides the QCP. 

Interpreting the phase diagram in these terms, one is 
immediately led to ask what the phase to the left of the 
QCP, i.e., for x < x c , might be. The theory would 
require a spontaneously broken symmetry beyond the 
global U(l) symmetry broken through superconductivity. 
In addition, as the fluctuations are assumed to determine 
the phase diagram up to temperatures of several hundred 
Kelvin, the characteristic energy scale of the correlations 
associated with this symmetry violation must be at least 
of the same order of magnitude. No definitive evidence of 
such a broken symmetry has been found up to now, even 
though several possibilities have been suggested. These 
include stripes 2 ^, a d-density wave^ 2 -, and a checkerboard 
charge density wave^ 3 .. 

The general consensus is that the low energy sec- 
tor of the three-band Hubbard model proposed for the 
CuO planes 2 ^ (see ^ below) reduces to a one-band 
t-t'-J model, with parameters t 0.44, t' w 0.06, 
and J ?s 0.128 (energies throughout this article are in 
e y-^ 25, 26,27, 28,29 _ However, two remarks are in order. 
First, the parameters are not exactly known, but can 
only be calculated approximately 3 ^. Second, for certain 
regimes of the phase diagram, CuO two-leg ladder studies 
have shown that the one-band and three-band description 
lead to qualitatively different result a 31 i 32 . 

For the undoped CuO planes, the formal valances are 
Cu 2+ and 2 ~. As the electron configuration of Cu 
atoms is [Ar] 3<i 10 4s 1 , this implies one hole per unit cell, 
which will predominantly occupy the 3d x 2_ y 2 orbital. As 
the on-site potential e p in the O 2p x and 2p y orbitals rel- 
ative to the Cu 3d x 2_ y 2 orbital is generally assumed to be 
of the order of e p = 3.6 (with = 0), and hence smaller 
than the on-site Coulomb repulsion Ua « 10.5 for a sec- 
ond hole in the 3d x 2_ y 2 orbital, additional holes doped 
into the planes will primarily reside on the Oxygens. The 
maximal gain in hybridization energy is achieved by plac- 
ing the additional hole in a combination of the surround- 
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FIG. 1: Orbital current pattern proposed by Varma. 

ing O 2p x and 2p y orbitals with the same symmetry as 
the original hole in the Cu 3d x 2_ y 2 orbital, which requires 
antisymmetry of the wave function in spin space, i.e., the 
two holes must form a singlet. This picture is strongly 
supported by data from NMR— and even more directly 
from spin-resolved photoemission2i In the effective one- 
band t-J model description of the CuO planes, these 
singlets constitute the "holes" moving in a background 
of spin 1/2 particles localized at the Cu sites. 

In contrast to this picture, Varma&iii has proposed 
that the additional holes doped in the CuO planes do 
not hybridize into Zhang- Rice singlets, but give rise to 
circular currents on O-Cu-0 triangles, which align into 
a planar pattern as shown in Fig. [TJ He assumes that 
the inter-atomic Coulomb potential Vpd is larger than 
both the hopping i p( j and the on-site potential e p of the 
O 2p orbitals relative to the Cu 2>d x 2_ y 2 orbitals, an as- 
sumption which is not consistent with the values gen- 
erally agreed on (see the list below ((T])). Making addi- 
tional assumptions, Varma has shown that the circular 
current patterns are stabilized in a mean field solution 
of the three-band Hubbard model. The orbital current 
patterns break time-reversal symmetry (T) and the dis- 
crete four-fold rotation symmetry on the lattice, but leave 
translational symmetry intact. The current pattern is as- 
sumed to disappear at a doping level of about x c ss 0.19. 
The phenomenology of CuO superconductors, including 
the pscudogap and the marginal Fermi liquid phase, are 
assumed to result from critical fluctuations around this 
QCP, as outlined above. 

Motivated by this proposal, several experimental 
groups have looked for signatures of orbital currents or 
T violation in CuO superconductors. While there is no 
agreement between different groups regarding the man- 
ifestation of T violation in ARPES studies^^, a re- 
cent neutron scattering experiment by Fauque et alr^ 
indicates magnetic order within the unit cells of the 
CuO planes. Their results appear to be consistent 
with Varma's proposal, and call the validity of the one- 
band models into question. In a recent article, Aji and 
Varma 3 -^ have mapped the four possible directions of the 
current patterns in each unit cell onto two Ising spins, 
and investigated the critical fluctuations. Within this 
framework, the couplings between and the transverse 
fields for these Ising spins effect whether or under which 
circumstances the model displays long-range order in the 
orbital currents. 
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We hence intended to estimate these couplings through 
numerical studies of finite clusters containing 8 unit cells, 
i.e., 8 Cu and 16 O sites, and periodic boundary condi- 
tions (which do not frustrate but should enhance the cor- 
relations). The total number of holes on our cluster was 
taken N = 9 (5 up-spins and 4 down-spins), correspond- 
ing to a hole doping of x = 1/8. We had hoped that the 
energy associated with a domain wall, which may be im- 
plemented through a twist in the boundary conditions, 
would provide information regarding the coupling align- 
ing the orbital currents in neighboring plaquets, while 
the splitting between the lowest energies for a finite sys- 
tem would provide an estimate for the transverse field. 
Together, this would account for a description of the sys- 
tem in terms of Ising-type variables where the quantum 
critical behavior could be analyzed. 

We find, however, that the current-current correlations 
in the ground state show no tendency to align the or- 
bital currents whatsoever. There is not even a context to 
speak of a coupling of these Ising variables — or, in other 
words, the couplings arc zero within the error bars of our 
numerical experiments. 



III. THREE-BAND MODEL FOR THE 
CUPRATES 



To begin with, we wish to study the three-band Hub- 
bard Hamiltonian H = H t + Htj with^ 

H t =X] e P n i^~*P d J2 ('C/' ; ..r ' /'U' / <, 

(id)," (i,j),<J,<? 



b) 



(1) 



where ( , } indicates that the sums extend over pairs of 
nearest neighbors, while di^ a and pj. a annihilate holes in 
Cu 3d x 2_ y 2 or O 2p orbitals, respectively. Hybertsen et 
al£L calculated t pd = 1.5, i pp = 0.65, U d = 10.5, U p = 4, 
Vpd = 1-2, and e p = 3.6, which is the first three-band 
model discussed below. 

In order to be able to diagonalizc (Q]) for a cluster of 
24 sites, i.e., 8 Cu and 16 O sites, with 5 up-spin and 4 
down-spin holes, we need to truncate the Hilbcrt space. 
A first step is to eliminate doubly occupied sites. This 
yields the effective three-band t-J Hamiltonian 

ff eff = P Q H t P G + Hj with 



H, = j pd J2( s^-sj - i ] + J DD V I Sf-S? - - 



where 
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pd \U d -e p T U p + e p J ' " pp U n ' 





FIG. 2: We define the phases q\ and qi acquired by trans- 
lations shown in a). For a system of eight unit cells, i.e., 8 
Cu and correspondingly 16 O, the phases can have the val- 
ues q\ = 7ii7r/2 and qi — ri2ir with n 1 £ {0,1,2,3} and 
n-2 £ {0, 1}. For the Brillouin zone depicted in b), qi and g2 
translate into x and y momenta by q x = qi and q y — q± — q% . 
Thus, the phase description (51,52) of the Y and M point is 
(0,0) and (7r,0), respectively. 

and the sums in Hj are limited to pairs where both neigh- 
bors are occupied by holes. If Pq only eliminates config- 
urations with more than one hole on a site, i.e., a pure 
Gutzwiller projection, the dimension of the S* ot = \ sub- 
sector is 164,745,504, which as such is above our capabil- 
ities. 

As a next step, we exploit the translational symme- 
tries on the cluster, a 4-fold symmetry in qi direction 
and a 2-fold symmetry in qi direction according to the 
conventions given in Fig. [2] (Throughout this article, 
we label the momenta by ((71,92) rather than (q x ,q y )-) 
The reduced dimension ~ 2 • 10 8 is still a considerably 
large Hilbert space. Thus, we take two further steps. 
First, we identify the ground state in the Brillouin zone. 
As it turns out, for all parameter choices discussed in 
this article, the respective ground state is either at the 
r or M point. Second, we apply the rotation symme- 
try, which commutes with the translational symmetries 
at and only at the M and V point. The system then 
becomes amenable to exact diagonalization. 

In order to identify the momentum of the ground state, 
we introduce two ways of truncating the Hilbert space: 
(a) We limit the maximal number of holes allowed in the 
O orbitals to iV™ ax . (b) We limit the maximal number 
of CuO links occupied with 2 holes to N^^. Truncation 
(a) serves as a good scheme when the on-site potential 
e p is large compared to other parameter scales, but is 
not practicable in other cases. Since truncation (b) pre- 
dominantly projects out states with high kinetic energy, 
we expect it to be insensitive to the value of e p . To 
check the validity of the truncations, we consider Hamil- 
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TABLE I: Ground state energies per unit cell for e p = 3.6, 
Vpd — 1.2 for the inequivalent points in the Brillouin zone and 
truncation (a). 
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FIG. 3: Current-current correlations (jk,k+xjl,m) multiplied 
by 10 2 and in units of £L for the ground state of ([2} with e p = 
3.6, V p d = 1.2 on a 24 site cluster (8 Cu = open circles, 16 O 
= filled circles) with PBCs. The reference link is indicated in 
the top and (due to the PBCs) bottom left corner. Except for 
the vertical lines, positive numbers indicate alignment with 
the pattern shown in Fig. [1] 



tonian (|2J with the parameter values by Hybertsen et al. 
listed above and calculate the ground state energies of 
the system (see Table HI and HT|) . 

Both truncation schemes yield similar results. The 
ground state is situated at the M point (tt, 0) of the 
Brillouin zone. At this point, it is possible to imple- 
ment a 4-fold rotation symmetry, which commutes with 
the translational symmetries. Thus, the dimension is re- 
duced by an additional factor of 4 and the exact state 
is accessible. We find the energies £^0,0) = —0.8513, 
£7(7r,o,7r/2) = —0.8570, and i^o,*-) = —0.8883, where the 
first two indices label the linear momenta (91,92), and 
the third labels the angular momentum under discrete 
rotations by 90°, q TQ t = n rot 7r/4, with n rot € {0, 1,2,3}. 
The ground state is hence in the (tt, 0, tt) subspace. 
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TABLE II: Ground state energies per unit cell for e p — 
3.6, Vpd = 1.2 for the inequivalent points in the Brillouin 
zone and truncation (b). 




FIG. 4: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of Q with e p = 3.6, V p d = 1.2 on 
a 24 site cluster (8 Cu = open circles, 16 O = filled circles) 
with PBCs. The Cu-O reference link is indicated in the top 
corner. 



IV. CURRENT-CURRENT CORRELATIONS 
AND MAGNETIC MOMENTS 

With the current operator for an 0-0 and a Cu-0 link 
given by 

jk.l = l -J^J2 (?l*Pk,<, - Pk,aPl,a) ( 3 ) 

(7 

and 

(7 

respectively, we evaluated the correlation function 
(jk.k+xji,m) with 0-0 and Cu-0 links as reference links 
for the exact ground state. The results are shown in 
Figs. [3] and [31 The correlations fall off rapidly and there 
is no indication of order. 

We now use the correlations to reach a conclusion re- 
garding the existence or non-existence of the orbital cur- 
rent pattern, indicated in Fig. [T] The numerical experi- 
ments for the finite cluster can, as a matter of principle, 
never rule out directly that a symmetry, in our case time 
reversal symmetry T, is violated. As only real parameters 
enter the Hamiltonian, the computed ground states are 
real by construction and do not allow for a direct indica- 
tion of time reversal symmetry breaking. If it were to ex- 
ist, the computed ground state would be a symmetric su- 
perposition of the different separately T violating ground 
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states, which itself is T symmetric again and described 
by a real wave function. The current-current correlation 
function, however, allows to put an upper bound on the 
size of the spontaneous currents: If a current pattern 
as sketched in Fig. [T] were to exist, the current-current 
correlations (jk,k+xjt,l+x) for links far away from each 
other in a rotationally invariant ground state should ap- 
proach i (x\jk,k+x\x) 2 , where \x) denotes a state with 
a spontaneous current pointing in x direction (the fac- 
tor i arises because by choosing our reference link in in- 
direction, we effectively project onto two of the four pos- 
sible directions for the current pattern) . From the values 
of 10 • (jk,k+x jl,l+x) for the four horizontally connected 
links in the center of Fig. El -0.0488, +0.0552, -0.0583, 
and +0.0552, which should all be positive if a current 
pattern were present, we estimate 10 2 • (jk.k+x) < 0.05 
and hence (x\jk,k+x \i) 2 < 1 • 10~ 3 as an upper bound for 
a current pattern we are unable to detect through the er- 
ror bars of our numerical experiment. (Throughout this 
article, currents are quoted in units of eV/h.) We now 
denote (x\j k M+x\x) by j pp . 

We roughly estimate the kinetic energy e pp per link 
associated with a spontaneous current j pp of this mag- 
nitude using j pp = n p v and £ pp = ^n p mv 2 with m = 
l/2t pp , where n p is the hole density on the Oxygen sites 
(rip = 0.14 for the state analyzed in Fig. [3j) , and obtain 
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A similar analysis with a Cu-0 reference link (shown 



in Fig. gl) yields with 10 2 



of 



(jk,k+x 
2 



< 1.0 and hence 



; < 10 2 (there is no factor i in this case) an estimate 



h 2 f 
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4£ p d v / n p"d 



< 5 ■ 10" 



(G) 



where we have determined via 8rid + 16n p = 9. 

Note that this energy of 3 (or 5) meV is not the conden- 
sation energy E c per unit cell, but a positive contribution 
to the energy of the current carrying state, which would 
have to be (more than) offset by other contributions (like 
the energy gain from aligning the circulating currents ac- 
cording to the pattern Varma proposed) if such a state 
were realized. We would expect the transition temper- 
ature T c of such a state to be of the order of the effec- 
tive coupling of the Ising spins introduced by Aji and 
Varma2£, while we would expect that E c < T c . 

Making contact to the quantities observed in experi- 
ment, we now derive an upper bound of the magnetic 
moment from the upper bound for the spontaneous cur- 
rents we obtained through numerics. The magnetization 
is related to the angular momentum of circulating elec- 
trons by M = —LiB,L z /h with L z = m e rv, where r de- 
notes the distance to the center of rotation and v the 
velocity of the electrons. 

With a Cu-Cu distance ao ~ 3.8A, the side lengths of 
the isosceles O-Cu-0 triangle are ^ for the two equal 




= eso/6\/2 



FIG. 5: Circulating current j pp on a neighboring O triangle 
according to the pattern in Fig. [pj Lengths are given in units 
of the Cu-Cu distance ao- 



legs and ^= for the third side, as shown in Fig. [5J The 
distance of the sides to the center-of-mass point are given 
by r p d = and r pp = j^jj=. Given the upper bound 

on the current-current correlations 10 2 {x\jk,k+x\x) 2 ^ 
0.1(^) 2 , we infer a bound for the particle current 



eV 

Jpp < 0.03 — . 



(7) 



Note that the units correspond to 1/time (leV/h = 
1.517 x 10 15 1/s). Assuming that each triangle is occu- 
pied by 1 hole only, this current corresponds classically to 
the inverse time the hole takes to go around the triangle 
once. 

Assuming further that the hole dwells equal amounts 
of time on each link defining the triangle, the velocity on 
the 0-0 link is given by 



ao i ( 1 1 
\/2 7 Ujpp 



3ao . 
V2 JPP 



and on the Cu-0 links by 



'J P d 



ao 



1 1 

3 jpp 



3a 



"Jpp- 



With L z 



wieT-pd^pd, we find an upper 



bound for the associated angular momentum 
£P p < 0.015 h. 



(8) 



As there are two current carrying triangles per unit cell 
in the pattern shown in Fig. [TJ we find an upper bound 
of the magnetic moment induced by the currents: 



Mceu = 2M A < 0.03 Mb- 



(9) 



The result is below the estimate of M w 0.05-0. 1/xb 
found experimentally by Fauque et al*£. (As we lower 
e p in Section IVTl below . however, the upper bound for j pp 
and jpd we are able to obtain from our numerical exper- 
iments increase, and becomes comparable to the experi- 
mental estimate.) 
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FIG. 6: An orbital current pattern proposed earlier by Varma, 
which has subsequently been ruled out by experiment. 

This provides an interesting perspective on the 
strength of the magnetic moment observed in the ex- 
periment. If a pattern as shown in Fig. [1] were respon- 
sible for this moment, the associated current would be 
in a range of ~ 0.05 — 0.1 eV/h, a signal strength which 
we do not observe in the numerics. Even if currents of 
this magnitude were present in the CuO layers, however, 
following the previous derivation, the energy associated 
with this current strength would only be in the range of 
10 — 30 mcV and thus too small to explain the phase dia- 
gram of the high-T c cuprates. This suggests that even if 
orbital current alignment were responsible for the mag- 
netic moment measured in experiments, these currents 
would not be sufficiently large to account for the strange 
metal phenomenology of high-T c superconductivity. 

The immediate conclusion we draw from our numerical 
results, however, is that if the CuO planes are adequately 
described by a three-band model with a set of couplings 
in the range we investigated, the antiferromagnetic order- 
ing observed by Fauque et alJ£ is not due to an orbital 
current pattern as shown in Fig. Q] We are hence led 
to ponder whether our models might be consistent with 
an alternative explanation of this experiment. The cor- 
relations we observe would be consistent with another 
orbital current pattern, shown in Fig. [SI which had been 
proposed by Varma earlier^Iii. This pattern, however, 
has been successfully ruled out by neutron scattering ex- 
periment^, and is not consistent with the experiment 
of Fauque et alr^. The correlations we observe, on the 
other hand, provide no indication that such a pattern is 
realized, as our system sizes are way to small to estab- 
lish the existence of any kind of long range order. We 
are merely not able to rule out order according to this 
pattern with our numerical data. Considering the possi- 
bility of this pattern being realized, however, does in any 
event not bring us closer to understanding the magnetic 
moment observed by Fauque et al^. 



SCALE = 100 



i 


'10.6941 } 

P -^2.721 


' 0.0198 

8 


* -0.0466 f 
/^HJ.695 


' 0.0198 ' 

2 


\ 




V- Hi 

0.0198 


-0.5919 } 

P-QWZ 
ir- HI 


V- Hi 

-0.1358 

7 


-0.0466 /* 
P -0.695 

i- HI 


2 


^ 


'-0.0466 } 

/ -0.69E 
|f— H| 


-0.1358 

2 


-0.0897 } 

/ -0.605 
— < 


-0.1358 

7 


► 


i 


0.0198 
*T7T^ 1 


-0.0466 

/^H).69£ 
HI 


-0.1358 

2 


-0.5919 

P -2.72C 


8 



10.6941 



0.019S 



-0.0466 



0.0198 



FIG. 7: Spin-spin correlations multiplied by 10 2 for the 
ground state of ((2j) with e p = 3.6, Vpd = 1.2 on a 24 site 
cluster (8 Cu = open circles, 16 O = filled circles). Due to 
PBCs the O reference site is in all corners of the plot and 
indicated by a big black circle. 

SCALE = 100 




FIG. 8: Spin-spin correlations multiplied by 10 2 for the 
ground state of ^ with e p = 3.6, Vpd = 1.2 with PBCs. 
The Cu reference site to the upper left is indicated by a big 
black circle. 



V. SPIN-SPIN CORRELATIONS 

If the observed magnetic moment is not due to orbital 
currents, what other possibilities are there? The only 
other explanation which comes to mind within the con- 
fines of our three-band models is antiferromagnetic order 
of the spins on the Oxygen lattice. The characteristic 



term driving the system into this kind of order would be 
the antiferromagnetic coupling J pp in our effective Hamil- 
tonian ((2J) . According to our intuition about the system, 
there is no reason to expect this kind of order, but as it 
is straightforward to obtain the spin-spin correlations for 
our finite cluster, we discuss them briefly in this section. 
We have evaluated the static spin-spin correlation of 
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spins at sites i and j given by 

Sij = (^SiSj^ , 

with an O or a Cu as reference site. The results are 
shown in Figs. [7] and respectively. For the O sublattice, 
we only observe a very small staggered spin correlation 
which falls off rapidly with distance on the scale of the 
lattice constant, regardless of whether we choose an O or 
a Cu as reference site. We do not observe any indication 
of antifcrromagnctic order on the O sites, as the required 
long range correlations are clearly absent. We will hence 
not discuss the spin-spin correlations any further and fo- 
cus on the current-current correlations instead. 

Regarding the experiment by Fauque et al^, we be- 
lieve that the explanation will require a model which goes 
beyond the projected three-band models studied here. 
The explanation might, for example, involve orbital cur- 
rents localized at the O atoms. 



VI. RESULTS FOR DIFFERENT PARAMETER 
SETTINGS 

We now consider various parameter choices for the 
three-band Hubbard model. The two parameters which 



We now turn to the different choices for the model pa- 
rameters. Except for e p and V* p d, we use the three-band 
Hubbard parameters calculated by Hybertsen et al£l. 
We label the different sections by the values for e p and 
V p d we specifically investigate. 



A. e p = 1.8, F pd = 1.2 

As we decrease e p from 3.6 to 1.8, the ground state 
switches from momentum (jr, 0) to (0,0), i.e.. from the 
M to the r point in the Brillouin zone (sec Tabic . 
Implementing the rotational symmetry yields -E(o,o,o) = 
-1.4551, Emp.Tr/2) = -1-4350, and E {0M = -1.4857 
(see Table HV|) . i.e., the ground state is at (0,0, tt). The 
current-current correlations for this state are depicted in 
Fig. [5] There is no evidence for a pattern along the lines 
of Fig. [TJ Instead, the correlations decrease rapidly with 
distance. The "fluctuations" or "noise" inherent in the 
finite size calculation are comparable to the preceding 
case e p = 3.6. The values of 10 2 • (jk,k+xji,i+x) for the 
four horizontally connected links in the center of Fig. [5J 
+0.0028, +0.0199, -0.0430, and +0.0199 indicate no rel- 
evant scale of correlations. Recall that the orbital current 
pattern shown in Fig. [TJ would require all the numbers to 
be positive. 



Varma^S, assumed to differ significantly from the values 
used above are the relative on-site potential e p , which he 
assumed to be small compared to all other energy scales, 
and the inter-atomic Coulomb interaction V~ p d, which he 
assumed to be the leading energy scale of the system 
generating the orbital currents. 

Thus, our strategy is as follows. Firstly, we succes- 
sively decrease e p from 3.6 to 1.8, 0.9, 0.4, and finally 0, 
and analyze each model as explained above. Secondly, 
for small values of e p , we double V p d from the standard 
value 1.2 to 2.4 and look whether this significantly in- 
fluences the system. In doing so, we implicitly sweep 
over a broad range of the charge transfer gap, which is 
dependent on V p ^ and e p( j. To begin with, we have com- 
puted the ground state energies per link for all incquiva- 
lcnt points in the Brillouin zone for all parameter choices 
we consider (see Table IIII[) . We found that the ground 
states are either situated in the T point or M point of the 
Brillouin zone. This is not surprising as it is plausible 
that the ground state does not carry any net momen- 
tum, and also consistent with the current pattern shown 
in Fig. [TJ For the states at the T and the M point, we 
then implement the rotation symmetry discussed above 
and compute the ground state energies per link in the 
respective subspaces exactly. The results are shown in 
Table [TVl 



r 





Vpd 


(0,0) 


(1,0) 


(*r,0) 


(o,0 


3.6 


1.2 


-0.8473 


-0.8356 


-0.8520 


-0.8359 


1.8 


1.2 


-1.4187 


-1.3888 


-1.4104 


-1.3925 


0.9 


1.2 


-1.7852 


-1.7508 


-1.7732 


-1.7569 


0.4 


1.2 


-2.0290 


-2.0011 


-2.0185 


-2.0066 


0.4 


2.4 


-1.6486 


-1.6381 


-1.6433 


-1.6387 


0.0 


1.2 


-2.2543 


-2.2427 


-2.2438 


-2.2407 


0.0 


2.4 


-1.9199 


-1.9265 


-1.9342 


-1.9261 



TABLE III: Ground state energies per unit cell for the dif- 
ferent momentum subspaces in the — 4 approximation. 
The numbers for the subspace with the lowest energies, cor- 
responding to the ground state of the full Hamiltonian, are 
underlined. 



B. e p = 0.9, Vp d = 1.2 

As we decrease e p further to 0.9, the lowest state re- 
mains in the T point, i.e., (0,0, tt). The ground state 
correlations are depicted in Fig. [TO] They fall off less 
rapidly with distance than for larger values of e p . Again, 
there is no indication of a current pattern. The values 
of 10 2 • {jk,k+xji,i+x} for the four horizontally connected 
links in the center of Fig. [TU1 -0.420, +0.0338, -0.0647, 
and +0.0338, indicate likewise no relevant scale of corre- 
lations. 
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v pd 


(0,0,0) 


(0A|) 


(0,0,7r) 


(7T,0,0) 


M,f) 


(7T,0,7r) 


3.6 


1.2 


-0.8544 


-0.8545 


-0.8843 


-0.8513 


-0.8570 


-0.8883 


1.8 


1.2 


-1.4551 


-1.4350 


-1.4857 


-1.4462 


-1.4448 


-1.4769 


0.9 


1.2 


-1.8440 


-1.8136 


-1.8672 


-1.8345 


-1.8304 


-1.8557 


0.4 


1.2 


-2.0986 


-2.0820 


-2.1149 


-2.0932 


-2.0859 


-2.1053 


0.4 


2.4 


-1.6864 


-1.6788 


-1.7007 


-1.6950 


-1.6741 


-1.6923 


0.0 


1.2 


-2.3269 


-2.3400 


-2.3389 


-2.3317 


-2.3280 


-2.3372 


0.0 


2.4 


-1.9508 


-1.9613 


-1.9647 


-1.9766 


-1.9543 


-1.9619 



TABLE IV: Exact ground state energies for the F point (0, 0, q TO t), and the M point (n, 0, q TO t), with additionally applied rotation 
symmetry (g rot = 3n/2 is degenerate to q Iot = 7r/2). The numbers for the subspace with the lowest energies, corresponding to 
the ground state of the full Hamiltonian, are underlined. 



+ convention for currents 



SCALE = 




+ convention for currents 



SCALE = 



FIG. 9: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of © with e p = 1.8, V pd = 1.2 
on a 24 site cluster with PBCs (reference link indicated in the 
top and bottom left corner). 

C. e p = 0.4, V pd = 1.2 

The correlations for the ground state, which remains 
at the (0,0, tt) point, are depicted in Fig. QT] We 
have now reached a parameter regime for which Varma 
proposed that a current pattern should occur: e p is 
small compared to Vpd, and V p d is of the order of the 
other scales. However, we still find no indication of 
a current pattern. The values 10 2 • (jk.k+xji,i+x) for 
the four horizontally connected links in the center of 
Fig. QH -0.0971, +0.0220, -0.0795, and +0.0220, remain 
small. Nonetheless, let us estimate the energy associ- 
ated with an upper bound for the currents as elaborated 
on above. We estimate 10 2 • (jk.k+x) 2 < 0.1 and hence 
(x\jk : k+x\x) 2 < 2 ■ 10 -3 as an upper bound for a uniform 




FIG. 10: Current-current correlations (jk,k+xjl,m) multiplied 
by 10 2 for the ground state of @ with e p = 6.9, V pd = 1.2 
on a 24 site cluster with PBCs. 



contribution according to a current pattern as depicted in 
Fig- HI As above, (£\jk,k+x\£) is denoted by j pp . For the 
kinetic energy e pp per link associated with a spontaneous 
current 



Jpp < 0-05 



eV 



we obtain with n r 



-pp 



0.24 

h 2 j 2 
Jpp 

4t pp Tip 



< 3 • 10" 3 . 



If a spontaneous currents were hence to exist, the energy 
associated with them would be too small to allow for an 
interpretation of the strange metal phase in terms of the 
quantum critical fluctuations around this phase. 
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+ convention for currents 



SCALE = 100 



+ convention for currents 



SCALE = 100 




FIG. 11: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of © with e p = 6.4, V pd = 1.2 
on a 24 site cluster with PBCs. 



D. e p = 0.4, Vpd = 2.4 

For e p = 0.4, we have also doubled the Coulomb 
repulsion between the Cu and O sites to V p d = 2.4, 
which then becomes the leading energy scale before the 
Copper-Oxygen hopping t p d = 1.5. The ground state 
remains at the T point of the Brillouin zone (see Ta- 
ble Compared to V p d = 1.2, the energies of the 
r and M ground states are now much closer to each 
other. The correlations for the ground state at (0, 0, 7r) 
are shown in Fig. [TJ] Again, there is no evidence of a 
current pattern. The values of 10 2 - (jk,k+xji,i+x) for the 
four horizontally connected links in the center of Fig. 1121 
-0.1501,-0.0487,-0.0399, and -0.0487, are compara- 
ble to the case Vpd = 1.2 discussed above. Note that the 
four horizontal links we consider are now aligned. Unfor- 
tunately, the sign of all four numbers is opposite to the 
sign required by the pattern shown in Fig. [TJ 



E. e p = 0, Vpd = 1.2 

This setting has already been discussed previously^. 
In a sense, e p = is the most favorable choice for Varma's 
mean field approach. It should be kept in mind, however, 
that e p must be positive and finite in the experimental 
system to account for the antiferromagnetic order in the 
undoped cupratcs. The ground state is now doubly de- 
generate and found at (0,0,±7r/2). There is no current 
pattern observable, but the values of 10 2 • (jk,k+xji.i+x) 




FIG. 12: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of (J2]) with e p = 0.4, V p d = 2.4 
on a 24 site cluster with PBCs. 

for the four horizontally connected links in the center 
of Fig. [T3 0.3726,-0.1483,-0.4123, and -0.1483, are 
larger than for any of the other parameter settings we 
investigated. The uniform current is still zero, but with 
larger "fluctuations" or "noise" due to the finite size of 
our system. We estimate an upper bound for a uni- 
form positive correlation 10 2 ■ (jk.k+x) 2 < 0.2 and hence 
(x\jk,k+x\x) 2 < 4-10~ 3 . As above, (x\j k ,k+x\x) is denoted 
by j P p, and the approximate kinetic energy e pp per link 
associated with a spontaneous current 

eV 

j pp < 0.06 — 

yields with ([5|) and n p = 0.34, 

£pp < 5 • 10- 3 . 

The correlations with Cu-0 reference link are shown in 
Fig. [TU A similar analysis yields with (jGj) 

e pd < 4- 10~ 3 . 

The upper bound of 5 meV we found for the energy asso- 
ciated with spontaneous currents corresponds to a tem- 
perature T curlGnt ~ 60K. To explain the strange metal 
phase, however, the energy scale responsible for the quan- 
tum critical fluctuations would have to extend at least to 
a range of several hundred Kelvin. 

Note that the upper bound for the magnetic moment 
per unit cell is now given by 

M ce u = 2M A < 0.06 mb, (10) 



9 




FIG. 13: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of © with e p = 6.0, V pd = 1.2 
on a 24 site cluster with PBCs. 

a value roughly comparable to the range of M w 
0.05-0. 1/xb found experimentally by Fauque et alH. 
This does not imply that these orbital currents exist, but 
only states that if we assume a three band model with 
e p = 0, we are not able to rule out an orbital current 
pattern as shown in Fig. [T] as an explanation for the ex- 
perimentally observed magnetic moment. On the other 
hand, even if the observed moments are due to such a 
pattern, the energies involved are too small to explain 
the phenomenology of the strange metal phase. 

F. e p = 0, Vpd = 2.4 

Finally, we double the repulsion V p d for e p = 0. The 
ground state is now at the M point, at {%, 0, 0) (see 
Tab. lIVp . The corresponding correlations are shown 
in Fig. 1151 Again, there is no evidence for a current 
pattern. The values of 10 2 • (jk,k+$jl,l+x) for the four 
horizontally connected links in the center of Fig. [TBI 
0.0909, -0.0885, 0.0865, and -0.0885, are rather small. 
We estimate an upper bound 10 2 • (jk,k+x) < 0.1 and 
hence {x\jk,k+x\x) 2 < 2 ■ 10~ 3 for the contribution from 
spontaneous currents. For the associated kinetic energy 
£ P p per link we find with n p = 0.36 

e pp < 2 • 10- 3 . 

Note that since Vpd is driving the orbital currents in the 
mean field theory proposed by Varma^, his theory would 
predict the model to show the strongest propensity to 




FIG. 14: Current-current correlations (jh,k+xjl,m} multiplied 
by 10 2 for the ground state of © with e p = 6.0, V p d = 1.2 
on a 24 site cluster with PBCs. The Cu-O reference link is 
indicated by a black arrow in the upper left corner. 

form current patterns for this choice of parameters. By 
contrast, we find no indication of such a propensity in 
our numerical experiments. 



VII. FINITE SIZE EFFECTS 

The numerical calculations we report here were per- 
formed on a cluster of 8 unit cells, i.e., 24 sites (8 Cop- 
per and 16 Oxygen). The question we wish to address in 
this section is whether any of the finite size effects might 
affect our overall conclusion. Due to the finite size, we 
have not been able to measure the current-current cor- 
relations at long distances, where they would accurately 
provide the size of a spontaneous current if such a cur- 
rent were to exist. Instead, we have only been able to 
establish upper bounds on such currents. These upper 
bounds, however, turned out to be small enough to allow 
us to rule out that an orbital current pattern as pro- 
posed by Varma and shown in Fig. [T] is responsible for 
the anomalous properties of the strange metal phase in 
the cuprates. 

There are, however, two other finite size effects which 
may limit the validity of our conclusion. The first is that 
the ground states in our systems have spin S = 1/2, as we 
have a total of 9 holes, corresponding to a doping of one 
hole away from half filling. The current carrying state 
proposed and investigated by Varma, by contrast, is in 
general a spin singlet. Could it be that the extra spin 
1/2 destroys the orbital current pattern in our numerical 
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+ convention for currents 



SCALE = 100 




FIG. 15: Current-current correlations (jk,k+xji,m) multiplied 
by 10 2 for the ground state of © with e p = 6.0, V pd = 2.4 
on a 24 site cluster with PBCs. 



experiments, while the currents would be present if we 
had an infinite system? We will argue now that any 
possible effect would not affect our conclusions. 

To begin with, assuming that no magnetic spin order 
is present, the ground state of the infinite system will 
be a spin singlet, regardless of whether the state car- 
ries a spontaneous current or not. This also holds for 
any finite system with an even number of electrons, even 
if long range antifcrromagnctic correlations in the spins 
were present. The ground state for a finite system with 
an odd number of electrons, as we have investigated in 
this work, will correspond to the ground state for an even 
number of electrons supplemented by an excitation which 
carries spin 1/2. This excitation will cost a finite amount 
of energy. The question relevant for the validity of our 
conclusion is whether the energy cost of this spin 1/2 ex- 
citation is higher for a current carrying state than it is 
for a state without currents, and if it is higher, by which 
amount. Since the term driving the spontaneous currents 
in Varma's analysis is Vpd, the energy associated with this 
term has to be lower in the current carrying state, while 
we would expect most other terms in the Hamiltonian, 
but in particular the antifcrromagnctic exchange terms 
J p d and J pp , to be slightly higher in energy. As the ex- 
citation energy for the extra spin 1/2 amounts mostly to 
an extra energy cost in J p( j and J pp , we expect that this 
excitation will cost less energy in the current carrying 
state than in the state without currents. In other words, 
unless the antifcrromagnctic exchange terms in |T]) were 
to contribute towards driving the system into a current 
carrying phase, the spin 1/2 excitation in our finite sys- 



tem would enhance the stability of this phase. The fact 
that we do not observe a current pattern in the presence 
of the extra spin 1/2 makes our conclusion even more 
robust than it would be without the excitation. 

To rule out any doubt completely, let us be unrea- 
sonable and assume that the antifcrromagnctic exchange 
terms do in fact enhance the systems propensity to de- 
velop spontaneous currents, and that the antifcrromag- 
netic exchange energy in the current carrying phase is 
maybe about 10% below the energy of the phase with- 
out the currents. For the parameter choice e p = and 
Vpd = 1.2, the antiferromagnetic exchange energy per 
unit cell is 

£J = i(* |if J |*o) = 122 mcV, 
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(11) 



with Hj given in To estimate the energy cost of the 
spin 1/2 excitation, we compare the ground state energy 
for 5=1/2 obtained for 5 up and 4 down spin holes 
with the ground state energy for 5 = 3/2 obtained for 6 
up and 3 down spin holes, and find an energy difference 
of A = 230 meV. Since the ground state energy for the 
different spin sectors should roughly be proportional to 
5 2 , the energy cost of the 5 = 1/2 excitation should be 
roughly l/8th of this difference, or about 30 mcV. The 
energy cost of the excitation per unit cell of our finite 
cluster is hence of the order of 4 meV. If we now assume 
that the energy cost for this 5=1/2 excitation increases 
by 10% if spontaneous currents are present (the cost in- 
creases because the excitation disturbs the current driv- 
ing antiferromagnetic correlations), the additional energy 
cost for the currents due to the extra spin 1/2 would be 
of order 0.4 meV per unit cell, or roughly 0.04 meV per 
link. If this energy cost were to destabilize the sponta- 
neous currents, the energy associated with them would 
be below the upper bounds estimated from the current- 
current correlations above. 

The second finite size effect we briefly wish to mention 
is that the special geometry of periodic boundary condi- 
tions might have an unintended influence on the system, 
maybe in that it stabilizes a state without currents which 
would not be stable in the infinite system. We have hence 
twisted and varied the boundary conditions in any way 
we could think of, but found that the correlations, and 
in particular the "fluctuation" or "noise level" due to the 
finite size which limits or ability to put upper bounds on 
the correlations, remained unchanged. 



VIII. CONCLUSION 

Let us summarize the results of our numerical studies 
on three-band Hubbard models for the cupratc planes in 
CuO superconductors. For the commonly accepted pa- 
rameter values as calculated by Hybertsen^I, we find no 
orbital current pattern as shown in Fig. [T] as well as no 
significant antiferromagnetic spin-spin correlation on the 
Oxygen sites. As we sweep over a considerable parame- 
ter regime of e p and Vpd and compute the current-current 
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correlations for the respective ground states, we likewise 
do not observe any evidence for a pattern as advocated 
by Varma's mean field analysis. Instead, we find that 
the correlations change quantitatively, but not qualita- 
tively, as we move through the parameter space. We also 
observe that there is no clear dependence of the correla- 
tions on Vpd, which is not consistent with what one would 
expect from Varma's analysis. We conclude that while 
we cannot rule out that orbital current patterns exist, 
we can rule out that they are responsible for the proper- 
ties of the strange metal phase or the anomalous normal 
state properties extending up to temperatures of several 
hundred Kelvin, as the energy associated with the spon- 
taneous loop currents would be less than 5 meV per link 
if such currents were to exist. We have assumed that the 
CuO planes arc adequately described by the three-band 
Hubbard model (fTJ), but we have allowed e p to be much 
smaller and V p d larger than generally agreed upon, and 



based our estimate for the upper bound of 5 meV on the 
for our purposes most unfavorable case e p = 0. 
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